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ABSTRACT 

The  feasibility  of  measuring  thermal  neutron  absorption  cross  sec- 
tion from  the  results  of  first-order  reactor  perturbation  theory  was 
investigated.   Techniques  were  developed  for  this  measurement  utilizing 
the  AGN-201  reactor  of  the  U.  S.  Naval  Postgraduate  School. 

These  techniques  were  evaluated  by  measurement  of  the  thermal  ab- 
sorption cross  section  of  indium  using  gold-197  as  a  standard.   The 
result  obtained  was 

194. 4;£  2.3  barns. 
A  comparison  with  the  most  recently  published  value  of  this  cross  sec- 
tion, 196ZI  5  barns,   demonstrates  the  feasibility  of  the  method  for 
obtaining  improved  precision  for  the  measurement  of  the  cross  sections 
of  selected  materials. 
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1.   INTRODUCTION. 

Thermal  neutron  absorption  cross  sections  now  reported  in  the  litera- 
ture have,  in  general,  precisions  exceeding  five  per  cent  although  measure- 
ments on  a  few  isotopes  have  been  made  to  much  better  precisions. 

One  method  of  measuring  neutron  absorption  cross  sections  is  suggested 
by  the  results  of  first-order  reactor  perturbation  theory  \_l].      This  method 
requires  the  comparison  of  changes  in  the  reactivity  of  a  reactor  caused 
by  the  presence  of  small  amounts  of  absorbing  materials  of  different  kinds. 

The  AGN-201  nuclear  reactor  of  the  U.  S.  Naval  Postgraduate  School 
is  designed  for  research  and  educational  purposes.   Because  of  the  small 
physical  size  of  its  core,  the  reactor  is  quite  sensitive  to  perturbations 
in  the  core  configuration.   The  purpose  of  this  investigation  was  to 
determine  the  feasibility  of  using  the  AGN-201  as  a  precision  device  for 
the  measurement  of  thermal  neutron  absorption  cross  section,  within  the 

framework  of  first-order  two-group  perturbation  theory. 

197 
A  foil  of  Au    was  chosen  as  a  standard  absorbing  material  because 

of  its  large  and  precisely  known  cross  section  (98.8t0. 3  barns),  the  con- 
venient half- life  of  its  reaction  product,  its  relative  ease  in  handling, 
and  its  natural  occurrence  as  one  stable  isotope.   A  foil  of  natural 
indium  was  chosen  as  an  "unknown"  for  the  purpose  of  investigating  this 
technique  because  of  its  large  cross  section,  and  the  relatively  simple 
decay  scheme  of  the  reaction  products. 

This  report  presents  a  brief  discussion  of  first— order  two-group 
perturbation  theory,  description  of  experimental  techniques  developed, 
results  and  error  analysis,  and  a  discussion  and  evaluation  of  the  methods 
employed. 


2.    THEORY. 


First  order  perturbation  theory  yields  an  expression  for  ^—^^P   ,  the 
change  in  reactivity  of  a  reactor  for  a  small  perturbation: 

„„     S0IV0JV 
^  =  J0/0dV 


where  ( /\    is  a  vector  representing  the  perturbed  reactor  flux,  (/)  '      is  a 


:or  flux,  (y\ 


1      i: 


vector  representing  the  unperturbed  adjoint  flux,  and  P  is  the  perturba- 
tion matrix.   The  integral  in  the  denominator  is  carried  out  over  the 
volume  of  the  reactor  as  a  whole  and  that  in  the  numerator  over  the  vol- 
ume of  the  perturbation.   The  number  of  energy  groups  determines  the 
number  of  components  of  both  vectors. 

For  the  two-group  case  the  vectors  may  be  written  as: 


0!<02  01 


and 


where  the  subscripts  1  and  2  identify  the  fast  and  slow  neutron  groups, 
respectively. 

For  this  case  and  when  the  perturbation  consists  only  of  small 
changes  in  the  absorption  cross  sections  in  a  small  volume  of  the  reactor 
caused  by  the  addition  of  a  poison, 


where   ^^~     =   tne  change  in  the  fast  neutron  average  macroscopic 

absorption  cross  section  and 

C<  -   =   the  change  in  the  slow  neutron  average  macroscopic 
absorption  cross  section, 

v  and  v9   =   the  average  speeds  of  the  neutron  groups,  respectively 
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Writing    A<?  :  Pf.     —«   ^o     ts   Che  —"vtty  with  a  void  in  the 

SMUvoJ,   1»    SI  a,  »d   Sfazb—   *.|    "*    *»*' 

scopic  cross  sections  of  the  poison  itself,  and  the  expression  for  A^> 

is: 

Uir^t^ti  ±  ^01  £~&W    (2) 

^°  "  Scat?,  +^/,0,)cJv 

If  the  effects  of  two  poisons,  identified  by  s  and  u,  are  being  considered, 
the  changes  they  produce  in  reactivity  may  be  compared: 

(p_p)      [m£!^  &  +^^^>tf  U^-(i> 

providing  the  denominator  of  equation  (2)  is  essentially  the  same  in  both 
cases.  Since  the  perturbation  is  assumed  to  extend  over  a  very  small 
volume,  r/~f'S  ~  i-u  /t(* 

or;5  (p0\     -    u, 


v/  011  =  v»tt  0 


<9  2- 


over  the  region  of  the  perturbation,  and  equation  (3)  may  be  more  con- 
veniently  written  as: 


t  (4) 


This  form  lends  itself  to  ready  interpretation.   The  integrals  represent 
the  slow  and  fast  absorption  rates  of  the  poisons.   The  quantity  in  curly 
brackets  is  a  weighting  function  which  expresses  the  relative  "importance" 
or  "worth"  of  fast  to  slow  neutrons  in  the  void  with  respect  to  overall 
reactor  performance. 

Assuming  the  poisons  are  homogeneous  and  defining 


y  =  t- 


equation  (4)  reduces  to; 


-?0)  _  zL0iv5+  ys.t  0,sv 


«      7yuv/u.     .      w  «-  ^    "wUi/U.  (5) 


(?u-R»;       s&ffv^+yzt.^v 


where  the  fluxes  are  now  volume -averaged.   Inspection  of  equation  (5) 
shows  that  it  relates  a  change  in  a  measurable  macroscopic  effect  (reac- 
tivity) to  certain  measurable  microscopic  processes  (thermal  and  epither- 
mal  neutron  absorption  rates).   Moreover,  if  poison  "s"  is  a  material  of 
known  cross  section  and  "u"  is  a  material  of  unknown  cross  section,  equa- 

tion  (5)  suggests  that  an  analytical  expression  for  <T     may  be  ob- 

/—  a  2- 

tained  using  ^\,_  as  a  standard,  provided  that  Y  and  the  necessary 

^  a  z. 

absorption  rates  and  reactivity  differences  can  be  determined. 

Appendix  I  develops  from  equation  (5)  expressions  for  the  ratio  of 
the  fast  neutron  absorption  rates  in  the  two  poisons  and  for  Y  in  terms 
of  experimentally  observable  quantities.   Substitution  of  these  back  into 
equation  (5)  thus  yields  the  following  formula  for  <£.         (equation  (A- 2)  of 
Appendix  I)  ; 


1"  =  i * 


where  the  previously  undefined  terms  are: 

_  Cd 

£-)     s   the  reactivity  with  an  empty  cadmium  cup  in  the  void 

region, 


Cd 
i       cadmium  cup,  and 


0  z      the  reactivity  with  the  foil  of  material  "i"  in  the 


i  for 

X  =  a  correction  factor/  the  foil  of  material  "i"  to 

account  for  fast  neutron  absorption  in  the  cadmium. 


3.    DESCRIPTION  OF  EXPERIMENTAL  TECHNIQUES. 

Inspection  of  equation  (6)  reveals  that  there  are  two  basic  types 
of  measurements  to  be  made;  reactivity  differences  for  various  configura- 
tions and  the  ratio  of  the  volume -averaged  thermal  fluxes  in  the  standard 
and  unknown.   This  section  describes  both  the  techniques  developed  for 
these  measurements  and  the  influence  of  geometry,  temperature  and  instru- 
mentation on  them. 

a)  Reactivity  Differences. 

There  are  two  basic  methods  of  obtaining  reactivity  differences 
with  the  AGN-201  reactor. 

One  is  from  the  measurement  of  the  change  in  the  stable  period 
caused  by  the  perturbation.   The  second  is  from  the  change  in  the  criti- 
cal position  of  the  control  rods. 

Preliminary  experiments  indicated  that  stable  period  measurements 
were  the  most  reliable,  because  the  reactor  operator  found  it  difficult 
to  determine  when  the  reactor  was  exactly  critical.   Even  though  the 
initial  timing  for  stable  period  measurement  was  accomplished  using  a  stop 
watch  and  observing  the  console  linear  micro -micro-ammeter,  the  results 
obtained  by  period  measurements  displayed  greater  reproducibility. 

Having  decided  to  use  stable  period  measurements,  it  was  necessary 
to  develop  reliable  instrumentation  for  timing  the  rise  in  reactor  power. 
Fig.  1  is  a  block  diagram  of  the  system  devised  for  measuring  times.  In 
order  to  obtain  a  larger  current  and  hence  less  fluctuation  than  In  the  out- 
put to  the  console  linear  micro-micro- ammeter,  a  Bfa  ion  chamber  was 
placed  in  an  access  port  of  the  reactor.   In  this  position  the  output 
current  was  approximately  100  times  greater  than  that  provided  by  the 
reactor  chamber.   The  signal  from  the  ion  chamber  was  fed  to  a  linear 
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micro-micro-ammeter  identical  to  that  of  the  reactor  console.   The  out- 
put of  the  micro-micro- ammeter  was  fed  to  a  positive-stop  micro- ammeter 
with  adjustable  stop.   When  during  measurement  the  micro-ammeter  read- 
ing first  corresponded  to  the  stop  setting,  the  battery  supplied  67.5 
volts  through  the  manual  switch  to  start  the  decade  scaler, which  counted 
10s000  cps  input  from  the  crystal- control- frequency  standard  meter. 
After  the  scaler  was  started,  the  scale  setting  on  the  micro-micro-ammeter 
was  increased  by  a  factor  of  ten  and  the  manual  switch  moved  to  the  "stop" 
position.   The  positive-stop  meter  then  stopped  the  scaler  when  the  power 
had  increased  by  a  factor  of  ten.   This  circuit  was  considered  precise 
to  about  5  milliseconds,  the  limiting  component  being  the  start-stop 
relay  of  the  scaler  itself,  and  this  did  not  introduce  any  significant 
uncertainties  into  the  decade  time  measurement.   The  conversion  from  the 
measured  decade  time  to  reactor  period  is  accomplished  by  dividing  dec- 
ade time  by  In  10.   The  table  of  the  inhour  equation  for  the  AGN-201 
reactor  of  Gans  (iff/  may  be  entered  with  the  stable  period  to  extract  the 
reactivity. 

The  following  techniques  were  used  to  measure  the  decade  times  for 
the  six  reactivities  appearing  in  equation  (6): 

(1)  The  control  rods  were  set  to  hold  the  reactor  at  or  near  criti- 
cal at  approximately  10  milliwatts  with  the  glory  hole  empty. 

(2)  A  sample  rod  containing  the  sample  holder  (discussed  later  in 
this  section)  was  rapidly  inserted  into  the  south  end  of  the  glory  hole 
(see  fig.  2).   Since  the  rod  loading  is  largely  polyethylene,  the  re- 
actor would  be  supercritical  for  all  perturbation  configurations  in  the 
void  region  of  the  sample  holder. 

(3)  The  power  was  allowed  to  rise  by  about  a  factor  of  thirty  before 
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commencing  decade  time  measurement  in  order  that  essentially  all  initial 
transients  would  have  disappeared. 

(4)  After  the  decade  time  measurement  had  been  completed,  the  sam- 
ple rod  was  removed,  and  a  special  cadmium  loaded  "killer"  rod  was 
promptly  inserted  to  return  the  reactor  to  10  milliwatts. 

(5)  The  killer  rod  was  removed  and  a  waiting  period  of  six  minutes 
observed  before  commencing  the  next  measurement  in  order  that  transients 
would  die  out.   The  entire  cycle  took  about  fifteen  minutes  and  the 
power  ranges  observed  were  from  about  10  milliwatts  to  five  watts.   At 
no  time  were  the  control  rod  settings  or  any  of  the  other  reactor  com- 
ponents changed.   In  any  sequence  of  measurements  the  configurations 
corresponding  to  \n  in  equation  (6)  were  done  first  to  determine  the 
most  reactive  configuration  for  safety  considerations.   Typical  periods 
ranged  from  about  fifteen  seconds  to  sixty  seconds. 

In  order  to  insure  that  the  only  changes  in  reactor  configuration 
from  run  to  run  occurred  in  the  sample  region,  a  special  sample  holder 
was  devised  and  a  clamp  was  used  on  the  sample  rod  to  insure  that  it  was 
identically  repositioned  in  the  glory  hole.   The  details  of  the  sample 

holder  are  shown  in  fig.  2.   The  holder  was  10  inches  long  and  consisted 

« 

of  two  mating  polyethylene  sleeves  with  inside  diameter  0.525  inch  in 
the  region  of  the  sample  and  0.500  inch  elsewhere.   Both  ends  were  in- 
ternally threaded.   Ten  one-inch  polyethylene  plugs  filled  the  bore. 
The  two  end  ones  were  threaded  so  that  they  might  be  used  to  butt  the 
plugs  tightly  against  each  other.   Thus  the  materials  in  the  sample  re- 
gion and  flux  monitoring  foils  at  other  positions  in  the  holder  could 
be  held  firmly  and  their  positions  reproduced  by  making  changes  in  only 
one  end  of  the  assembly. 
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With  half- inch  diameter  foils,  preliminary  runs  made  it  apparent 
that  a  10-mil  thickness  would  be  required  to  provide  sufficiently  large 
reactivity  changes  for  precise  determinations  of  the  differences  in 
reactivities  with  and  without  the  sample  present. 

The  holder  in  the  \       condition  contained  all  materials  which  were 
required  to  be  in  the  holder  during  flux  determinations  except  the 
sample  itself.  (The  specifications  of  these  materials  will  be  discussed 
in  Section  b).   Briefly  they  are:   (1)  at  the  sample  position,  which  is 
at  the  flux  peak  one  inch  north  of  core  center,  an  incident-flux  monitor 
foil  surrounded  by  a  gold  ring,  and  (2)  at  the  flux  peak  one  inch  south 
of  core  center,  a  power- time  monitor.   The  r  s  and  P  measurements  were 
made  with  the  10-mil  gold  and  indium  foils,  respectively,  positioned 
within  the  gold  ring  and  held  snugly  against  the  incident-flux  monitor. 
The  cadmium  runs  were  made  with  the  same  loading  except  that  in  each 
case  the  gold  ring  with  its  contents  was  placed  inside  a  cadmium  cup 
20  mils  thick. 

Early  in  the  experiment^  on  reactivity  determination,  it  became 
apparent  that  changes  in  reactor  temperature  were  affecting  reproduci- 
bility of  the  determinations.   This  is  not  surprising  since  the  static 
temperature  coefficient  of  the  reactor  is  large  and  was  reported  by 
Harvey  and  Ferguson  [7 J  to  be  -0.02837.,/  C.   Therefore,  changes  of  just 
a  few  tenths  of  a  degree  could  be  very  significant  in  determining  reac- 
tivity of  a  given  sample  from  a  series  of  measurements.   An  attempt  was 
made  to  measure  small  temperature  changes  in  order  to  correct  the  reac- 
tivity measurement.   Since  the  tele thermometer  installed  on  the  reactor 
console  cannot  be  read  to  less  thanO. 3  C,  an  amplified  sensitivity  was 
required.   This  was  obtained  by  use  of  the  recorder  output  of  the 
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tele thermome ter.   By  use  of  the  circuit  shown  in  fig.  3,  small  changes 
were  recorded  on  the  graphic  recorder. 


1. 

34  V 

O- 
Yellow  Springs 

Instrument  Co. 

YSl  Model  425C 

Tele thermome ter 

G- 

10  K 

AAA 

t 

-o 

Varian  Associates 
Model  G-14 
Graphic  Recorder 

(0-1  millivolt 

scale) 

FIGURE  3 
CIRCUIT  FOR  MEASURING  SMALL  CHANGES  IN  TEMPERATURE 

This  recorder  had  a  full  scale  deflection  of  approximately  0.7°C  when  used 
in  this  circuit.   Although  it  consistently  indicated  small  temperature 
changes  in  the  reactor  during  the  reactivity  runs,  the  magnitudes  of  the 
changes  could  not  be  determined  with  enough  certainty  to  be  useful. 

The  form  of  the  reactivity  terms  in  equation  (6)  made  it  possible 
to  minimize  the  significance  of  temperature  changes.   The  six  reactivi- 
ties that  appear  fall  into  two  natural  groups  of  three  each: 

Cd   _  Cd   _  Cd 
(  \u'  \s'     Vo)  and  (Vu   »   Is   '   Vo  )'   The  three  members  of  each 

group  are  then  used  for  two  reactivity  differences.   Assuming  that  a  sig- 
nificant temperature  change  would  not  occur  within  the  time  required  to 
measure  a  full  group  ('>-  one  hour),  it  was  then  possible  to  obtain  reac- 
tivity differences  within  each  set  of  measurements  with  a  minimum  sys- 
tematic error  introduced  by  temperature  changes.   The  precision  desired 
in  the  experiment  could  then  be  obtained  by  repeating  the  sets  of 
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measurements  several  times.   In  order  to  further  minimize  the  effects  of 
temperature,  the  reactor  was  run  at  five  watts  for  about  twenty  minutes 
in  an  attempt  to  establish  temperature  equilibrium  before  beginning 
period  measurements.   In  addition, the  timing  sequence  during  successive 
runs  was  carefully  reproduced.   These  procedures  gave  consistent  results 
during  the  final  reactivity  difference  determinations. 

b)  Ratio  of  the  Volume -averaged  Thermal  Fluxes. 

The  factor^  /0   in  equation  (6)  may  be  obtained  from  the  ratio 
of  the  thermal  activity  of  the  standard  foil  to  the  volume-averaged 
thermal  flux  in  the  unknown  foil.   To  obtain  the  former,  standard  cadmium 
ratio  techniques  were  utilized  and  the  activities  corrected  for  differences 
in  radiation  time  and  reactor  power  level  by  using  the  activities  of  a 
gold  power-time  monitor  located  elsewhere  in  the  reactorj^J   was  obtained 
by  relating  the  thermal  activation  of  a  thin  gold  "incident  flux"  monitor, 
adjacent  to  the  unknown,  to  the  spatial  distribution  of  the  thermal  flux 
in  the  unknown,  measured  in  a  separate  experiment.   Activities  were  ob- 
tained by  relative  gamma  counting  and  in  all  cases  were  extrapolated  to 
time  of  removal  from  the  reactor,  corrected  for  dead  time,  background, 
and  when  necessary,  self-shielding.   The  gold  counting  data  were  reduced 
using  the  program  of  Appendix  II. 

All  counting  was  conducted  using  a  3"  X  3"  Harshaw  Nal(Tl)  scintil- 
lation detector,  a  type  12512  Dumont  photomultiplier  tube  with  a  cathode 
follower,  a  Hamner  non-overload  amplifier  (model  N-300),  a  Hamner  high- 
voltage  power  supply  (model  N-401),  and  a  Hamner  Decade  Scaler  (model 
N-240).   Gamma  counts  were  taken  in  the  integral  mode  using   low-voltage 
selection  to  reduce  the  effects  of  gain  shift.   An  aluminum  cap,  3"  I.D. 

and  1/4"  thick,  was  fitted  to  the  scintillation  detector  head.   A  small 
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"well,"  1/32  inch  deep  by  0.502  inch  in  diameter,  was  machined  into  the 
center  of  the  surface  of  the  cap;  this  device  served  a  twofold  purpose. 
The  aluminum  acted  as  a  shield  to  avoid  possible  interference  from  beta 
emissions.   The  "well"  served  as  a  positive  positioning  aid  for  the 
foils  counted,  a  necessary  step  since  the  detection  efficiency  is  quite 
sensitive  to  change  in  solid  angle.   It  was  also  necessary  to  place  a 
1/4-lb  lead  weight  atop  the  foils  while  counting, because  even  slight 
curvature  in  an  activated  foil  caused  a  difference  in  the  counting  rates 
obtained  by  turning  the  foil  over,  of  the  order  of  3-4  standard  devia- 
tions. 

All  foils  were  counted  to  a  statistical  precision  of  0.1%  or  better 
on  each  side.   Corrections  were  made  for  background,  which  in  general 
was  not  negligible,  and  dead  time  loss.   (The  dead  time  for  the  counting 
equipment  was  determined  by  the  double-pulse  method,  checked  by  an  oscil- 
loscope and  found  to  be  1.75  £  .05  J-l   sec.) 

Because  many  of  the  experiments  involved  counting  many  foils  over 

a  total  of  several  hours,  precautions  were  taken  to  monitor  the  stabil- 

137 
ity  of  the  counting  equipment.   A  C_   source,  mounted  on  a  copper  plan- 

chette  and  with  a  count  rate  comparable  to  that  of  the  activated  foils, 
was  counted  before  and  after  each  experiment.   No  data  were  used  for 
which  the  Cs  count  rates  differed  by  more  than  two  standard  deviations. 
The  counting  equipment  was  preconditioned  before  measurements  were  made 
by  counting  an  activated  foil  for  several  minutes. 

Since  relative  ^-counting  was  employed,  no  corrections  were  made 
for  counter  efficiency,  backing  gain,  or  geometrical  factors.   A  correc- 
tion for  a   -      self-shielding  was  incorporated  where  necessary;  this 
correction  is  of  the  forra£  5j  :   £  ~\~p^ 


-■*4f>* 
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where  jV(    z     self  shielding  attenuation  coefficient  =  0.19  cm  /gm 

P  =   density  of  the  material  being  counted  =  19.296  gm/cnr  for 
rolled  gold 

t      =   thickness  of  the  material  (cm) 

Although  the  correction  is  small  for  the  0.5  mil  gold  foils  employed, 

it  is  quite  significant  for  the  10-mil  foils.   The  ratio  of  these  cor- 


rections, which  will  be  utilized  in  the  calculation  for  (/)      /7j     ,   is 
1.045  ±  0.001. 


0SM 


To  relate  (/)        to  the  activities  of  the  bare  and  cadmium-covered 


0 

standard  foil  irradiated  at  different  times,  first  write,  for  the  bare 
case,  the  average  count  rate  corrected  to  time  of  removal  from  reactor. 


where  Ns  =   the  average  observed  count  rate  =   A  Ns/(l-e   c) 

N   =   the  number  of  counts  due  to  gold  in  time  t 

r\      =  decay  constant  of  Au  198  =  2.9  7462  X  10"6/sec. 

ts  r   the  time  from  removal  of  the  foil  from  the  reactor  to  the 
w 

beginning  of  the  counting  interval 

tr  =   the  time  of  irradiation 

£  z    /■  self  shielding  correction 

i\    z     overall  counting  efficiency  excluding  self  shielding 

L   a  a  dimensionless  parameter  accounting  for  the  actual  power 
level  of  the  irradiation. 


Similarly,  the  average  count  rate,  corrected  to  the  time  of  removal  from 
the  reactor,  for  the  gold  power- time  monitor  during  the  same  irradiation 
is 
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■s     — 


t-f,e'*--[£MY<,®,V'Ii-e"'Il£] 


-  [£,  f  ♦£,  0;  j  =  iftf  k.  0;  *  £  en 

S        P 

where  M  and  M  are  the  masses  of  the  standard  and  time-monitor  foils, 
respectively. 

Similarly,  for  the  cadmium  shielded  irradiation 

For  these  experiments  the  standard  foil  was  enclosed  in  a  20-mil  thick 
cadmium  cup  which  was  carefully  sealed  to  prevent  thermal  neutron  pene- 
tration. 

Combining  these  results  gives 

®--(tMM„w--ah)[L2):*£,0: 

\/J 2     can  be  related  to  the  thermal  activation  of  a  thin  gold  inci- 
dent-flux monitor  adjacent  to  the  unknown  foil  in  the  following  way. 
Let  Ki  be  the  ratio  of  (7)  to  the  thermal  flux  at  the  interface  between 


the  unknown  foil  and  the  incident-flux  monitor, and  K0  be  the  ratio  of 

m 

2 


the  flux  at  the  interface  and^y   ,the  average  thermal  flux  in  the  monitor 
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itself.   Then 


(/J       can  be  determined  from  counting  experiments  in  the  same  manner  as 
was  j£y2  •   Expressed  in  terms  of  activities ,  j£y   is 


where  the  prime  distinguishes  the  power- time  monitor  for  these  two  irradia- 
tions from  that  for  the  standard  foil. 

Finally,  since  the  two  power- time  monitors  were  essentially  identi- 
cal, 

The  gold  foils  used  were  0.500"  diameter.   In  addition  to  the  10- 
mil  gold  standard  three  other  thicknesses  were  employed;  2-mil  gold  foils 
served  as  "power- time"  monitors,  0.5-mil  foils  were  used  as  "incident- 
flux"  monitors,  and  1-mil  foils  were  used  in  a  "sandwich"  technique  for 
determining  Ko  as  described  below.   Foils  were  cut  using  a  light  duty 
metal  punch.   Copper  backing  to  minimize  foil  damage  was  utilized.   Foil 
masses  were  determined  by  weighing  on  a  Mettler  type  H-16  electronic 
balance  and  are  precise  to  £  0.03  mg.   Prior  to  weighing,  all  foils  were 
cleaned  using  a  50/50  dry-ether/ethanol  solution. 
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It  was  felt  desirable  to  have  the  sample  at  one  of  the  two  flux 
"peaks"  (fig.  2)  and  the  power-time  monitor  at  the  other.   In  these  rela- 
tively "flat"  regions  of  the  flux  profile,  small  imprecisions  in  reposi- 
tioning of  the  sample  rod  would  not  cause  an  appreciable  error  in  the 
relative  flux  "seen"  by  either  monitors  or  sample. 

To  ascertain  that  the  cadmium  perturbation  of  the  ambient  flux 
would  not  significantly  affect  the  power-time  monitor,  the  following 
technique  was  adopted: 

Two-mil  gold  foils  were  inserted  between  the  blocks  of  the  poly- 
ethylene holder  assembly  described  earlier,  one  foil  in  each  of  six 
positions,   when  the  entire  holder  assembly  and  sample  rod  were  fully 
inserted  into  the  "glory  hole"  the  foil  positions  corresponded  to  one, 
two,  three  and  four  inches  south  of  core  center,  and  three  and  four  inches 
north  of  core  center.   The  sample  position  was  at  one  inch  north  of  core 
center.   Two  irradiations  were  made,  the  first  being  conducted  with  the 
cadmium  cup  in  the  sample  position,  and  the  second  with  the  sample  posi- 
tion empty.   Construction  of  the  sample  holder  guaranteed  precise 
repositioning  of  the  foils  for  the  two  runs. 

It  was  expected  that  the  ratio  of  specific  activities  of  foils  at 
corresponding  positions  would  be  identical  within  statistics   in  the 
regions  unaffected  by  the  cadmium,  while  those  ratios  for  foils  in  the 
neighborhood  of  the  cadmium  would  show  a  departure  from  the  mean  ratio. 
The  ratio  for  the  foils  at  the  desired  position  of  the  monitor  (i.e. 
one  inch  south  of  core  center)  showed  no  such  departure,  and  hence  was 
selected  as  an  appropriate  position  for  the  power  time  monitor. 


An  interesting  effect  was  noted  during  early  attempts  to  determine 
these  ratios.   When  the  specific  total  activity  ratio  of  two  foils  of 
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In  order  to  determine  K  and  K«,  a  "sandwich"  technique  was  util- 
ized.  A  stack  of  thin  foils (of  gold  for  K2)  was  irradiated  with  the 
incident- flux  monitor  at  one  face.   The  specific  thermal  activation  of 
each  foil  was  determined,  and  plotted  as  a  function  of  position  in  the 
stack.   (The  data  for  indium  are  shown  in  fig.  4  which  appears  in  Section 
41 

A  gold  ring,  5  mils  thick,  1/16"  wide  and  slightly  less  than  0.5" 
in  diameter  was  constructed  (see  Section  5  for  details),  and  into  this 
ring  the  "sandwich"  was  inserted.   The  purpose  of  the  gold  ring  was  two- 
fold:  (1)  to  eliminate  over -activation  of  the  inner  foils  caused  by 
small  separations  between  their  edges,  and  (2)  to  provide  some  physical 
support  for  the  sandwich  while  encased  in  the  cadmium  cup  during  the 
fast  neutron  activation  runs. 

The  gold  sandwich  consisted  of,  in  order,  the  incident-flux  monitor, 
two  additional  half-mil  gold  foils,  and  nine  one-mil  gold  foils.   The 
indium  sandwich  consisted  of  the  gold  incident-flux  monitor  and  five  two- 
mil  indium  foils. 

Ki  was  obtained  from  the  graph  of  fig.  4  and  the  computed  total 
specific  activities  of  indium.  oiT>)( incident)  ,  which  is  directly  propor- 
tional to  the  thermal  flux  at  the  interface  between  the  indium  foils  and 
the  incident  flux  monitor,  was  obtained  by  extrapolation  of  the  data 
points. 


considerably  different  masses  was  compared  with  that  of  two  foils  of 
nearly  the  same  mass,  the  heavier  foil  was  "colder"  than  the  lighter  foil, 
In  the  cases  noted,  the  two  foils  differed  in  mass  by  about  107».   The 
apparent  difference  in  specific  activity  was  approximately  17c.   This 
anomaly  is  probably  due  to  the  difference  in  behavior  of  the  epithermal 
and  thermal  flux  through  the  volume  of  the  foil.   An  approximate  calcu- 
lation of  the  relative  flux  depression  and  their  effect  on  specific 
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Ko  was  obtained  from  a  similar  plot  of  the  gold  data. 
The  gold  ring  and  the  incident  flux  monitor  thus  become  an  integral 
part  of  the  holder  assembly  for  all  reactivity  and  flux  measurements. 


activity  through  two  gold  foils,  differing  in  thickness  by  107o,  using 
the  known  thermal  cross-section  and  effective  resonance  integral  of 
gold,  substantiates  this  explanation.   Therefore,  the  foils  used  in  all 
subsequent  traverse  experiments  and  as  power- time  monitors  had  masses 
matched  to  within  0.25  mg. 
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4.    RESULTS  AND  ERROR  ANALYSIS. 

a)  Reactivity  Measurements. 

The  results  of  the  reactivity  measurements  are  expressed  first  in 
terms  of  the  four  reactivity  differences  of  equation  (6);  these  are  then 
combined  with  the  cadmium  correction  factors  to  yield  the  overall  result: 

0^  CJ    Cd, 


Six  measurements  were  made  of  the  set  (  p  ,  p  ,  0   )  and  three  of  the 

Cd    Cd    Cd  u        o 

set  (  P   ,   P  ,  P  ).   The  means  of  the  four  reactivity  differences 

and  their  variances  were  calculated  assuming  the  differences  were  norm- 
ally distributed.  It  is  noted  that  by  treating  each  reactivity  differ- 
ence as  a  random  variable,  no  allowance  was  made  for  the  fact  that  I   and 

Pcd 

\    each  appear  in  two  terms.   This  introduces  a  correlation  in  the  es- 
timated variances.   Because  of  the  form  of  D,  the  estimated  variances 
will  be  larger  than  the  true  ones. 

The  values  of  the  cadmium  correction  factors  are  taken  from  refer- 
ence [_3J.   Their  standard  deviations  are  estimated.   Column  1  of  Table  1 
lists  these  results  and  their  estimated  standard  deviations.   Reactivity 
differences  are   expressed  in  per  cent.   Column  2  shows  the  contri- 
bution of  each  term  to  the  variance  of  D  which  was  calculated  using 
Gaussian  error  propagation. 
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TABLE  I 
RESULTS  OF  REACTIVITY  MEASUREMENTS 


Column   1 

C< 
( 

slumn   2 

6 
X  10      ) 

Pu" 

Po 

0.02571±  0.00012 

47 

Ps- 

Po 

0.020 18  ±0.00014 

88 

c 

-     ff 

0.00395+0.00007 

19 

esd 

-     tf 

0.00267  ±10.00002 

1 

rs 

1.015    +.005 

3 

*u 

1.095    +.005 
1.224  ±  .013 

5 

0 

163 
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b)  Ratio  of  the  Volume -averaged  Thermal  Fluxes. 

Difficulty  was  experienced  in  obtaining  reproducible  results  for 
the  sandwich  experiments;  this  is  discussed  at  length  in  Section  4.   The 
results  from  the  most  reliable  data  were  used  to  calculate  K^  and  K~. 
Fig.  4  is  a  plot  of  the  data  for  Ki.   The  estimated  precision  in  the 
specific  activities,  which  included  counting  statistics  and  the  uncer- 
tainties due  to  mass  measurements  only,  were  about  0.10%.   The  estimated 
precisions  of  Ki  and  Ko  include  uncertainties  due  to  the  graphical  analy- 
ses used  and  to  other  random  effects  intrinsic  to  the  sandwich  technique. 
The  values  are: 

KL  =  0.9890  £  0.0054 

K2  r  0.9816  ±  0.0014 
The  counting  rates  for  the  standard  foil  and  the  incident  flux  moni- 
tor are  precise  to  approximately  0.17».   Foil  mass  measurements  were  known 
to  about  0.03  mg.   Combining  the  uncertainties  in  all  factors  gives 


0: 


s 

?■    -    1.0610±  0.0060 
u 


The  precision  is  dominated  by  the  uncertainties  in  K,  and  K2« 

c)  The  2,200  m/sec  Microscopic  Cross  Section  of  Indium. 

Equation  (6)  may  be  rewritten  in  terms  of  the  microscopic  neutron 
absorption  cross  sections  of  the  standard  and  unknown  materials  averaged 
over  the  thermal  neutron  spectrum  of  the  AGN-201  reactor  at  20.0  C. 
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where  W  and  W  are  the  atomic  weights  of  the  standard  and  unknown  ma- 

s      u  ° 

terials,  respectively,  and  T  is  the  absolute  temperature  of  the  reactor 
during  the  experiments.  The  results  for  the  factors  in  this  expression 
and  the  estimates  of  their  standard  deviations  are: 

CT&2  98.8  ±  0.3  barns      [j*~} 

Mg  614.34  ±  0.05 

^  234.16  ±  0.05  mg 


w 

wu 


Ws 


&'/&" 


D 
T 

Combining  these  quantities  gives 


114.82 

C8J 

197.0 

M 

1.0610  £  .0060 

1.224  ±   .013 

294.5  ±  0.3  °K 

u 
a2 


Q*   ;   =   194.4  t    2.3  barns. 


Assuming  that  the  energy  dependences  of  the  gold  and  indium  cross  sec- 

__  u 
tions  in  the  thermal  range  are  not  significantly  different,  Q"*  ?  may 

be  compared  to  196"X  5,  the  published  0.0253-ev  value  of  the  cross  sec- 
tion for  indium  (j^U*  * 
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5.    DISCUSSION  AND  CONCLUSIONS. 

The  result  obtained  for  the  cross  section  of  indium  is  of  consider- 
ably better  precision  than  that  reported  in  current  literature  providing 
no  hidden  systematic  error  is  present.   It  is  apparent,  however,  that 
even  more  precise  determinations  are  possible  with  further  refinement 
of  these  techniques.   The  two  measurements  that  limit  the  precision  of 
the  determination  in  this  experiment  are  bare  reactivity  differences  and 
the  flux  depression  factor  K1K0  which  relates  the  average  flux  in  the 
incident  flux  monitor  to  the  average  flux  through  the  sample  foil. 

The  determination  of  reactivities  may  prove  quite  difficult  to 
improve.   Although  the  AGN-201  is  well  suited  to  measurements  of  reac- 
tivity changes  caused  by  samples  of  small  physical  size,  it  is  also  very 
sensitive  to  temperature  changes  which  may  be  caused  by  a  change  in  en- 
vironmental conditions  as  well  as  internal  changes.   A  considerable  num- 
ber of  trials  were  conducted  in  order  to  determine  how  to  achieve  maxi- 
mum temperature  stability  and  how  to  measure  small  temperature  changes. 
Although  a  technique  was  developed  to  minimize  the  effects  of  temperature 
change,  no  positive  means  of  measuring  these  changes  was  achieved. 

One  method  of  further  minimizing  the  temperature  effect  would  be 
to  take  sets  of  measurements  consisting  of  only  the  two  reactivities 

occurring  in  each  of  the  reactivity  differences.   This  would  require 

_         Cd 
determination  of  £sQ   and  (O     each  time  a  difference  in  which  they  are 

used  is  measured. 

The  method  which,  if  accomplished,  would  provide  much  better  con- 
trol over  the  experiment  would  be  the  design  of  a  device  capable  of 
precise  detection  of  small  temperature  changes  and  a  careful  measurement 
of  the  dynamic  temperature  coefficient  of  reactivity  at  low  powers.   The 
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sensor  for  such  a  device  must  be  investigated,  however,  as  apparent  tem- 
perature changes  of  small  magnitude  as  observed  on  the  specially  designed 
temperature  recorder  detected  by  the  sensor  of  the  tele thermometer  were 
not  accompanied  by  the  immediate  reactivity  effects  expected  with  this 
change;  usually  the  reactivity  effects  were  observed  at  some  later  time. 
Furthermore,  use  of  the  static  temperature  coefficient  of  Harvey  and 
Ferguson  did  not  give  consistent  correlations. 

A  possible  source  of  systematic  error  was  present  in  the  reactivity 
determinations  because  the  foils  were  subjected  to  bending  and  warpage 
when  being  fitted  inside  the  gold  ring.   The  fact  that  the  sample  was 
changed  for  each  successive  run  increased  the  opportunity  for  such  dam- 
age to  occur.   Foil  handling  could  be  reduced  to  a  minimum  if  proper  tem- 
perature corrections  could  be  made,  thus  allowing  repeated  measurements 
of  the  same  configuration. 

As  with  reactivity  measurements,  the  sandwich  technique  used  for 
determination  of  average  thermal  flux  in  the  unknown  could  be  consider- 
ably improved. 

Several  early  experiments  were  devoted  to  application  of  this  method 
to  a  gold  sandwich  of  the  same  size,  but  composed  of  nine  l~mil  gold 
foils  in  the  sandwich  body  with  two  half-mil  foils  at  the  surface  adja- 
cent to  the  half-mil  monitor.   In  all  instances  save  one,  however,  the 
technique  failed;  the  cadmium- shielded  irradiations  produced  anomalously 
high  activities  among  the  lighter  foils,  leading  directly  to  an  apparent 
thermal  depression  at  these  positions.   This  was  probably  caused  by  some 
separation  of  these  small  foils  from  the  body  of  the  sandwich  during 
the  cadmium  irradiation.   The  single  success  involved  a  sandwich  com- 
posed of  four  two-mil,  one  one-mil,  and  three  half -mil  foils  (this  series 
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is  that  from  which  K2  was  computed) . 

To  fully  appreciate  the  origins  of  this  condition,  it  is  necessary 
to  review  the  details  of  preparing  the  sandwich-ring  assembly.   The 
sandwich  is  contained  by  a  5-mil  X  1/16"  gold  ring,  of  slightly  smaller 
diameter  than  that  of  the  sandwich  itself.   This  ring  is  fabricated  in 
the  following  manner:   a  strip  of  gold  about  1/16"  wide  is  cut  from  a 
sheet  of  5-mil  stock  in  the  shape  of  an  elongated  cross,  the  arms  of 
which  have  a  span  of  about  3/16".   The  head  of  the  cross  is  beaten  down 
to  a  thin  taper  and  trimmed  to  the  width  of  the  cross  body.   The  strip 
is  then  rolled  around  a  polyethylene  cylinder  which  has  previously  been 
cut  down  at  the  end  by  about  two  mils  from  an  original  .500"  diameter. 
The  body  of  the  cross  is  pulled  back  tightly  across  the  tapered  head, 
the  arms  fold  up  and  over  and  are  then  pressed  down,  locking  the  assembly 
in  the  shape  of  a  ring.   Extra  material  is  trimmed  away.   The  ring  may 
be  then  carefully  worked  up  the  cylinder  until  it  is  held  by  only  about 
1/4  its  thickness.   Foils  may  be  inserted  at  this  time,  pressing  them 
down  with  undercut  polyethylene  blocks  so  that  the  edges  of  the  foils 
bind  tightly  against  the  sides  of  the  ring.   If  accomplished  with  great 
care,  the  whole  assembly  may  be  lifted  free  of  the  matrix  and  inserted 
in  its  cadmium  shield.   The  same  loading  procedure  is  followed  after  the 
cadmium  irradiation  is  completed  and  counted. 

This  cadmium  irradiation  is  the  most  difficult  phase.   While  within 
the  cadmium  cup,  the  foils  have  no  positive  agent  pressing  them  tightly 
together.   Instead,  the  foils  are  supported  only  by  the  pressure  of  their 
edges  against  the  sides  of  the  ring.   It  is  apparent  that  error  in  load- 
ing technique  might  allow  a  poor  purchase,  especially  among  the  thin 
foils,  allowing  the  crucial  thinner  foils  to  work  loose  during  irradiation, 
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"Hype ractivat ion"  would  result.   The  fact  that  this  effect  was  not  noted 
on  the  single  successful  run  may  be  attributed  to  the  fact  that  two-mil 
foils  were  utilized  in  the  body  of  the  sandwich  during  this  experiment, 
rather  than  one  mil.   These  heavier  foils  may  have  allowed  a  more  uni- 
form loading,  and  the  additional  stiffness  behind  the  outer  monitors  may 
have  prevented  their  working  loose  during  irradiation. 

The  indium  sandwiches  suffered  from  their  own  problems  peculiar  to 
the  nature  of  the  material.   Because  indium  tends  to  be  somewhat  adher- 
ent, no  trouble  was  experienced  because  of  separation  of  the  foils  in 
the  sandwich.   However,  because  indium  is  also  extremely  soft,  foil 
damage  was  most  difficult  to  avoid,  especially  in  the  form  of  weight 
loss  (at  the  edges)  during  loading.   This  necessitated  reweighing  the 
foils  after  each  irradiation. 

Edge  damage,  not  only  to  indium,  but  to  gold  also,  may  have  contri- 
buted to  the  poor  results;  it  is  almost  certain  that  some  foils  were 
somewhat  marred  about  the  circumference  during  one  of  the  loadings. 

The  fact  that  the  gold  ring  could  conceivably  slip  relative  to  the 
sandwich  proper,  allowing  one  face  to  be  more  in  the  neutrons  "shadow" 
of  the  protruding  gold  cylinder  than  the  other,  might  have  introduced 
error.   Pains  were  taken  to  prevent  such  an  occurrence  by  machine  tool- 
ing the  blocks  and  by  careful  visual  inspection  of  the  sandwich  before 
and  after  each  irradiation. 

In  summary,  it  is  felt  that  this  method  of  flux  depression  measure- 
ment is  feasible  in  principle,  but  that  a  significant  improvement  in 
both  precision  and  reliability  could  be  achieved  by  the  development  of  a 
better  mechanical  scheme  for  loading  and  holding  the  foils  in  the  sand- 
wich. 
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In  conclusion,  the  authors  believe  that  the  feasibility  of  the 
methods  developed  for  cross  section  measurements  has  been  demonstrated 
and  that,  although  the  reported  results  reflect  a  considerable  improve- 
ment on  other  methods  of  thermal  cross  section  measurement  for  certain 
isotopes,  the  precision  of  1.2%  attained  in  this  investigation  can  be 
bettered  with  a  minimum  of  further  experimentation. 
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APPENDIX  I 
EXPRESSION  OF  EQUATION  (6)  IN  TERMS  OF  EXPERIMENTALLY  OBSERVABLE  PARAMETERS 
Suitable  algebraic  manipulation  permits  (6)  to  be  solved  explicitly 

tor£. 


u 


(A-l) 


Besides  the  reactivity  differences,  <£?u  -f^Q)   and  ^s_P0)>  an<*  tne  ratio  of 
the  volume -averaged  thermal  fluxes  in  the  two  materials,  the  quantities 
which  must  be  measured  are  the  ratio  of  the  fast  neutron  absorption  rates 
in  the  unknown  and  standard  materials,  the  ratio  of  the  fast  to  slow  neutron 
absorption  rates  in  the  standard  material,  and  the  relative  worth  of  fast 
and  slow  neutrons  (Y) . 


£a;^svs 


The  ratio  of  the  fast  neutron  absorption  rates  can  be  obtained 
by  using  cadmium  cups  to  filter  out  the  thermal  neutrons.   If  the 
void  reactivity  of  equation  (5)  refers  to  a  void  inside  a  cadmium  cup, 
equation  (5)  becomes: 

U  — 3U      u  v  ^s 


f  r-?r .  i 


^-j  Cd      s-\  Cd  ^-    sbs      s  k/u 

Ks  "To        Z«i0i  v  0 


iT 


? 


Where: 
Cd 


u 

Cd      . 
s 


Reactivity  observed  for  cadmium- shielded  unknown  material 
Reactivity  observed  for  cadmium- shielded  standard  material 
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£2  =  Reactivity  observed  for  the  empty  cadmium  shielding. 

z:w  V-      b 

Then  the  ratio  <—- may  be 

.  s  —7  s  s 

£..  pi  v 

yu    /pCd     pCd 
written  as  Q_     /  \   u        \  o 

Cd  ^  Cd 


rie."-^ 


r 


Substituting  the  cadmium-covered  standard  material  for  the  unknown 
in  equation  (5)  : 

rj  .  Cd  --/Cd  Cd      _  Cd^Cd*  Cd     .  srns   s   1 

fs  f .       £*  ^ v  +  Y1,  <p,  v  +Y1,0, v  T 

Similarly,  substituting  the  empty  cadmium  cup  for  the  unknown  in  equa- 
tion (5)1: 


The  validity  of  this  expression  depends  upon  the  approximate  equality 


or 


of  the  quantity   J  ((^o(  $/  "J"  (3  T  02.)dV  witn  tne  cadmi-um  CUP  *n 

out  of  the  reactor.   A  rough  calculation  utilizing  the  known  spatial  de- 
pendence of  the  flux  for  the  two  cases  indicated  that  the  percentage  dif- 
ference is  about  0.13%  which  is  smaller  than  other  errors  in  the  determina- 
tion of  Y. 
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Cd 


ri  cd  sA  <-- 

Note  that  (/)      differs  from  UJ  because  the  epithermal  flux 

in  the  cadmium  is  affected  by  the  presence  of  the  standard  material. 

Inverting  these  two  expressions,  subtracting  and  solving  for 
Y  gives : 


£ 


5 


a-7- 


t  0 


I  /-f)  i  H£rfV*-£ 


•;  .  5 1 


xr^ 


t  &  v? 


Q~ 


Rewriting  the  third  term  in  the  denominator  as 


a 


and  approximating 


tf£  vci 


I 


0 


Cd 


*Y 


and       ^i    0i      V  by      ]f    -  J 


the  term  may  be  approximated 


r 


as 


r 


which  for  gold  is  about  0.00022.   Therefore,  the 


third  term  of  the  denominator  may  be  neglected  in  comparison  with 
the  others  within  the  experimental  precision  expected. 
Within  these  approximations  Y  may  be  expressed  as: 
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U.^U  v  ,U 


Substitution  of  these  expressions  for/ ,    _ — ^— ]    and  Y  into 


2i,yt  Vs 


equation  (A-l)  gives: 


U 


±    -  2 


■Q.3. 


.5 


t  V5!  fCfi,-?a)-ytf-?f j)  >  ca-2) 


5 


Tj  (JQ-W^WW) 


U.  \/UL 

2 


It  is  to  be  noted  that  the  ratio  of  the  fast  to  slow  neutron  absorption 
rates  does  not  appear  in  this  equation. 
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APPENDIX  II 
COMPUTER  PROGRAM  FOR  REDUCTION  OF  GOLD  ACTIVATION  DATA 

The  following  Fortran  60  computer  program,  written  by  Associate 
Professor  E.  A.  Milne,  was  utilized  for  computation  of  the  gold  foil 
specific  activities. 


DIMENSION  D(100),C(100),R(100) 
5  FORMAT  (1H1,/21X,6HNUMBER,11X,7HWAITING,7X8HC0UNTING6X9HC0RRECTED, 
18X7HRATE  AT/3X7HFOIL  ID11X9HOF  COUNTS8X4HTIME10X4HTIME10X4HRATE13X 
26HTW  =0/) 

10  FORMAT  (3E10.5) 

11  FORMAT  (A8,A7,F10.0,3F10.5) 

12  FORMAT  (A8,A7,F15.0,4F15.5) 

13  FORMAT  (/20X16HTRAVERSE  NUMBER  12/) 
I  =  1 

J  =  0 

M  ■  1 

WRITE  OUTPUT  TAPE  3  ,  5 

READ  INPUT  TAPE  2,  10,  FLAM,  TAU,  BG 

15  READ  INPUT  TAPE  2,  11,  A,B,D(I) ,TC ,TW,W 
IF  (D(I))  99,99,16 

16  C(I)  =  D(I)/TC 

C(I)  =  C(I)/(1.0  -  C(I)*TAU)  -  BG 

R(I)  =  C(I)*EXPF(FLAM*TW)/W 

IF  (I-M)  20,18,20 
18  M  ■  M   12 

J  =  J   1 

WRITE  OUTPUT  TAPE  3,  13,  J 
20  WRITE  OUTPUT  TAPE  3,  12,  A,B,D(I) ,TW,TC,C(I) ,R(I) 

1  =  1   1 

GO  TO  15 
99  STOP 
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